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Abstract 

We consider the neutral inclusion problem in three dimensions which is to prove if a 
coated structure consisting of a core and a shell is neutral to all uniform fields, then the 
core and the shell must be concentric balls if the matrix is isotropic and confocal ellipsoids 
if the matrix is anisotropic. We first derive an over-determined boundary value problem 
in the shell of the neutral inclusion, and then prove in the isotropic case that if the over¬ 
determined problem admits a solution, then the core and the shell must be concentric balls. 
As a consequence it is proved that the structure is neutral to all uniform fields if and only if 
it consists of concentric balls provided that the coefficient of the core is larger than that of 
the shell. 
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1 Introduction 

The purpose of this paper is to prove that the coated inclusions neutral to all uniform fields in 
the isotropic medium are concentric balls in three dimensions. The coated inclusion is depicted 
by ( D , Q.) where D and O are bounded domains with Lipschitz boundaries in M d (d = 2,3) such 
that D C fh Here, D represents the core and Q\D the shell. The conductivity (or the dielectric 
constant) is a c in the core and a s in the shell (er c y cr s ). If the structure (D,Q) is inserted into 
the free space with conductivity a m where there is a uniform field — V(a • x) = —a for some 
constant vector a, then the field is perturbed in general. But for certain inclusions the field is 
not perturbed, in other words, the field does not recognize the existence of the inclusion. For 
example, the coated inclusion is made of concentric balls with specially chosen conductivities 
(confocal ellipsoids if a m is anisotropic), one can see the uniform field is not perturbed. The 
inclusion with this property is called a neutral inclusion (or neutrally coated inclusion) and the 
neutral inclusion problem is to show that the inclusions of concentric balls (or confocal ellipsoids) 
are the only coated inclusions neutral to all uniform fields. 
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Nos. 2010-0017532 and 2013R1A1A1A05009699, and by Grant-in-Aid for Scientific Research (B) (jf 26287020) of 
Japan Society for the Promotion of Science. 
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Let a denote the conductivity distribution of the medium so that 


o- c 

in D, 

0's 

in 12 \ D, 

Om 

in R d \fl 


( 1 . 1 ) 


Here we assume that a c and a s are constants (or isotropic matrices), but a m is allowed to be 
anisotropic symmetric matrix. We consider the following problem: 

f V • a Vu = 0 in R d , 

{ ’ ( 1 . 2 ) 

I m(x) — a • x = 0(|x| “) as |x| —>• oo, 

where a is a constant vector. The term u(x) — a • x depicts the perturbation of the potential 
due to insertion of the coated inclusion (D,Q). If the potential is not perturbed, namely, 

u(x) — a ■ x = 0 in \ H, (1.3) 


the coated inclusion ( D , fi) is said to be neutral to the field a. If ( D , H) is neutral to the field 
ej for j = 1,... ,d, where e 3 is the standard basis of R d , then (H,H) is neutral to all uniform 
fields. 

Much interest in neutrally coated inclusions was aroused by the work of Hashin and Shtrik- 
man [7] and Hashin [Bj. They showed that since insertion of neutral inclusions does not perturb 
the outside uniform field, the effective conductivity of the assemblage filled with coated inclu¬ 
sions of many different scales is a m . We refer to m for developments on neutral inclusions in 
relation to the theory of composites. Another interest in neutral inclusions has been aroused 
in relation to invisibility cloaking. The neutral inclusion is invisible from the probe by uniform 
fields as observed in m Recently, the idea of neutrally coated inclusions has been extended 
to construct multi-coated circular structures which are neutral not only to uniform fields but 
also to fields of higher order up to N for a given integer N (2]. It was proved there that the 
multi-coated structure combined with a transformation dramatically enhances the near cloaking 
of [13]. Cloaking by transformation optics was proposed in [IT) (and [5j). 

As mentioned before, concentric balls (or disks) are neutral to all uniform fields by choosing 
<j c , u s and a m properly (a m is isotropic). Confocal ellipsoids (or ellipses) are also neutral to all 
uniform fields if a m is anisotropic P2J (see also section 0. Then a question arises naturally: 
are there any other shapes which are neutral to all uniform fields? In two dimensions there are 
no other shapes: if a coated inclusion (Z),fl) is neutral to all uniform fields in two dimensions, 
then D and H are concentric disks (confocal ellipses if a m is anisotropic). This is proved when 
<7 C = 0 or oo in m and when a c is finite in [TO] . In this paper we consider the neutral inclusion 
problem in three dimensions. We emphasize that the methods in [ TOl m use powerful tools 
from complex analysis such as conformal mappings and harmonic conjugates, which cannot be 
applied to three dimensions. It is worth mentioning that there are many different shapes of 
coated inclusions neutral to a single uniform held as shown in two dimensions in ©EH- 

We first show that if ( D , H) is neutral to all uniform fields in three dimensions and if cr c > cr s , 
then the following problem admits a solution: 

Aw = k in H \ D, 

< Vw = 0 on <9f2, (1.4) 

Vw(x) = Ax + d on dD , 

where k{^ 0) is a constant, A is a symmetric matrix, and d is a constant vector. We emphasize 
that this is an over-determined problem because Vw is prescribed on the boundaries. The 
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problem, which is of independent interest, is to prove that if (11.41) admits a solution in three 
dimensions, then D and 0 are confocal ellipsoids. If D and 0 are confocal ellipsoids, then (11.41) 
admits a solution and A should be either positive or negative-definite depending on the sign of k 
(see section El). So a part of the problem is to show that A is either positive or negative-definite. 
In two dimensions it is proved in EQ3 that if (11.41) admits a solution then D and Q are confocal 
ellipses (concentric disks if A is isotropic). However, the proof there is based on the powerful 
result that there is a conformal mapping from Ll\D onto an annulus. So it cannot be extended 
to three dimensions. The condition er c > cr s , which is not natural, is required because of a 
technical reason for the derivation of ([1.41) in subsection 12.21 Even though we do not know how 
to do so, it is likely that the condition can be removed. 

In this paper we solve the problem partially as the following theorem shows. 

Theorem 1.1 Let D and 0 be bounded domains with Lipschitz boundaries in R 3 with D C f2. 
Suppose that Ll\D is connected. vm admits a solution for A = cl for some constant c where 
I is the identity matrix in three dimensions, then D and O are concentric balls. 

As a consequence, we obtain the following theorem. 

Theorem 1.2 Let D and H be bounded domains with Lipschitz boundaries in R 3 with Dell. 
Suppose that dD is connected and R 3 \ D is simply connected. If a m is isotropic, <j c > a s and 
(D, f l) is neutral to all uniform fields, then D and H are concentric balls. 

This paper is organized as follows. In section El we show that if ( D , 0) is neutral to all 
uniform fields then (11.41) admits a solution. In section El we construct a solution to (ll.4|) when D 
and f l are confocal ellipsoids. Section E] is to prove Theorem ll.il In section El we formulate the 
problem (11.41) using Newtonian potentials and relate the problem with a known characterization 
of ellipsoids. 


2 Derivation of the over-determined problem 


In this section we derive (11.41) out of the neutral inclusion problem. We will do so only in three 
dimensions since (|1.4j) has been derived in two dimensions m ■ We assume that dD is connected 
and R 3 \ D is simply connected. 

Suppose, after diagonalization, that 

&m — diag[<7 mi i, <7 mi 2, cr m,3] • (2.1) 


Let Uj, j = 1, 2, 3, be the solution to 


f V • aVuj = 0 
\ Uj(x ) — Xj = 0(\x\~ 2 ) 



as \x\ —>• oo. 


( 2 . 2 ) 


The structure being neutral to all three fields means that Uj (x) — Xj 
Let 

Wj = l. Uj 

where 

J' = 1,2,3, 

(Js 

and w = (wi, W 2 ,ws) t (T for transpose). Set also 


0 in M 3 \ Q for j = 1, 2,3. 

(2.3) 


B = diag [l//3i, l//3 2 , !,%]■ 


(2.4) 


We will show the following: 
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(i) Vw is symmetric and divw is constant, and hence there is a function '</; in Q\D such that 

w = VV> and Atp = TrB + 1 in P \ D. (2-5) 

(ii) w(x) = cox + d, x € D for some constant Co and constant vector d (under the assumption 
that a c > a s ). 

We emphasize that it is in (ii) where the condition a c > a s is required. 

Once we have (i) and (ii), then we can show that m has a solution. In fact, since Uj = Xj 
on dO, we have 

V^(x) = Bx on dO. 

Note that V^x) = cox + d on dD. Now define 

rc(x) := ^(x) — ^x • Bx. (2.6) 

Then w satisfies (11.41) with k = 1 and A = col — B. We emphasize that if a m is isotropic, so are 
B and A. 


2.1 Proof of (i) 

Let us first deal with the case when 0 < a c < oo. Denote by v = ( n\,n2,n^) T the outward unit 
normal vector field to dO or dD. Note that the solution Uj (j = 1,2,3) to (|2.2D satisfies the 
following transmission conditions on two interfaces : 


and 


diLj 

— u j\- =0, 

diin 


- Uj \ — = 0, (T s 


dv 


- v. 


duj 

ls ~dJJ 

duj_ 

dv 


= 0 on dP 


= 0 on dD 


(2.7) 


( 2 . 8 ) 


where + denotes the limit from outside and — that from inside of P or D. If (D, P) is neutral 
to Xj, then Uj(x) — Xj = 0 in M 3 \ P, so we see from (12.71) that 


Uj | — — Xj) (Tg 


du j 

dv 


= CTmjrij on dP. 


(2.9) 


In other words, Uj is the solution to the following over-determined problem: 


V • aS/uj = 0 
duj 

u ,=x 3 , — 


a 


V s 


Tli 


in P, 
on dP. 


( 2 . 10 ) 


Let Vj £ C 2 (P). Then we see from the divergence theorem and (12.81) that 


/ 

Jd n 


duj_ 

dv 


dvj 

V >~ U ^ = - 


/ UjAvj + l 

Jfl\D J dD 


dUj_ 

dv 


Vn ~ U 


dvj 

‘ 3 lhJ 


ln\D 


UjAvj + 


JQ\D 

= — UjAVi + 


(j. 

Vs 


UjAvj + (— - l) f Vj+ f 

\Vs J J dD dv - Jg 

/ Vn 7 - • Vm — 

Jd Jl 


du i 
i 9 d dv 


Vi — u 


d^ 

3 dv 


UjAvj 


v. 


Vs 


- 1 


>D 


\7uj ■ Vfj. 
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On the other hand, we see from (12.91) that 


/ 

Jan 


du j 

dv 


dvj 

f &m, 

£ 

1 

e 

'■o. 

c»| 
^ l> 

II 

Jan cr s 


1 a m,j 


V cr s 


-rij v j Vj 


dvj_ 

dv 


- 1 


9vi 

dy 


n u yj Jn 


~ / Vj Av j- 


Equating two identities above we obtain 


[ {Vj - Uj)A Vj + a f V Uj • Vvj - Pj f = 0, j = 1, 2,3 (2.11) 

Jn Jd Jn a Uj 


for Vj £ C* 2 (9), where a and (3j are defined for ease of notation to be 


cr, 


a = -1 and /L = 


(7 


cr. 


m,] 

cr. 


- 1. 


( 2 . 12 ) 


Let Wj be defined by iVj := -j^Uj as in (12.31) . Then (12.lip can be rephrased as 


1 f f dv 

(l5-yj~ w j) Av j + a Vwj-Vvj- 7rf = 0, j = 1,2,3. (2.13) 

in Pj Jd Jn a Vj 


Summing (12.131) over j = 1, 2,3 we have 


8Vn 


I _ w j) A v 3 +«/ Vu h • Vv j - / Y1 dt = 0 

In Pj Jd “( Jn °yj 


j= 1 




for v-j £ C* 2 (9). If we use vector notation w = (wi,W 2 ,ws) t and v = {vi,V 2 ,vs) t (T for 


transpose), then the above identity can be rewritten as 

/ (By — w) ■ Av + a / Vw : Vv — / divv = 0. 
Jn J d Jn 


(2.14) 


Here and afterwards A : B denote the contraction of two matrices A and B, i.e., A : B = 
Y,a ij b lj = TV(A r B). 

Let r be the fundamental solution of the Laplace operator in R 3 , i.e., 


r ( x ) := “TT7’ x ^°- 

47r|x| 


(2.15) 


Let Vj( y) = T(x — y) for a fixed x £ 9. Since Avj (y) = <5(x — y), by applying the divergence 
theorem over 9 \ B e (x) for sufficiently small e (B e (x) is the ball of radius e centered at x) we 
see from (12.13|) that 

u>j(x) = ^-Xj + a [ Vwj(y) ■ V y r(x — y)dy + -^—Nq(x), x £ 9, j = 1,2,3, (2.16) 

Pj Jd oxj 

where Nq is the Newtonian potential on a domain 9, i.e., 

JVo(x) := [ F(x-y)dy, x £ R 3 . (2.17) 

Jn 


Let 


fj{x):= Vwj(y) • V y r(x — y)dy, j = 1,2,3, 


ID 


5 














and let f = (/i, /h, fz) T ■ Note that fj is harmonic in R 3 \ D, and (|2.16l) can be rewritten as 

w(x) = af(x) + V Qx • Bx + Afo(x)^ , x G 11, j = 1,2,3. (2.18) 

For any fixed x £ I 3 \ U, let 

^'( y ) = ^- r ( x - y ), 3 = 1 , 2 , 3 . 

Then divv(y) = —A y T(x — y) = 0 and Av(y) = 0 for y G 12. So we see from (|2.14l) that 


ID 


Vw : Vv = 0, 


(2.19) 


and hence 

divf(x) = f VVw ? (y)-V^-r(x-yWy= f Vw : Vv = 0, xGR 3 \12. 

Jd j dxj J D 

Since fj is harmonic in R 3 \ D, (12.191) holds for all x G R 3 \ D. 

Again fix x G R 3 \ 12. Let {i,j, k} be a permutation of {1, 2, 3} and let 

d d 

y) = — r (x-y), Vj(y) = r(x - y), v k = 0, y € 12. 

Then, Av = 0 and divv = 0 in 12. So we have from (12.141) 

J^ Vitfi(; y) • V^-T(x - y)dy - Vw 3 {y) • V^-T(x - y)dy = 0, 
which implies that 

difM) = fi ( x ) 

for all x £ R 3 \ 11 and hence for all x £ R 3 \ D. Moreover, since R 3 \ D is simply connected, by 
the Stokes theorem there is tp such that 


( 2 . 20 ) 


Because of (12.191) . we have 

Let 


f(x) = Vp(x), x£R 3 \D. 


A <p(x) = 0, x£R 3 \D. 


1 


i/j(x) = atp(x) + -x • Bx + Nq(x), x G 12 \ 12. 
Then, we have from (12.18|) and (12.211) 

w(x) = V^(x), x£ll\D, j = 1,2,3. 
Since A Nq(x) = 1 for x G 12, we have from (|2.22l) that 

Aijj(x) = Tr B + 1, xgH\D. 


( 2 . 21 ) 

( 2 . 22 ) 

(2.23) 

(2.24) 

(2.25) 


So far we have shown that Vw is symmetric, divw is constant, and (12.51) holds when a c is finite. 
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We now assume that a c = 0. In this case the problem (12.101) becomes 

A Uj = 0 in 12 \ D, 

^ = 0 on dD, 

dv 


(2.26) 


dllj &m,j ;w. 

Uj = Xj, —— =- -rij on dll. 


07 


So, we see in a way similar to ( 12 . 111 ) that 
/ Vj Av 3 - / 

J £1 J Q 


in\D 


uiAvj - L u ^ - p* L 


dvj_ 

n d Vj 


= 0 


(2.27) 


for all Vj e C 2 (fl). So we obtain a representation of the solution similar to (12.161) : 


w 


( x ) = -j -xj - wj(y)^ - y)dtr(y) + ^-Afo(x), xgH\D. 


(2.28) 


So, we infer in the exactly same way as in the previous sections that Vw is symmetric and divw 
is constant, and there is a function ip such that (12.51) holds. 

Suppose that a c = oo. In this case the problem (12.101) becomes 


A Uj = 0 


u 


j = 7 j (constant) 


in Q\D, 
on dD, 


9uj cr,n.} .... 

Uj = Xj, —— = - -rij on oil. 

ov cr. 


(2.29) 


The constant 7 j is determined by the condition 

f d u j 
J dD 9v 

We then obtain similarly to (12.111) 


= 0. 


[ vA v i ~ [ u A v i + [ 7^7 v i ~ 7 j [ A '7 - Pj 

JQ JD\D J dD JD 


'n dVj 


= 0 


(2.30) 


for all Vj e C 2 (fl). We then obtain a representation of the solution similar to (12.161) : 


w 


1 C r)in ■ r) 

(x) = — Xj + — 1 (y)r(x-y)do-(y) + — JV n (x), x <E 11 \ D. (2.31) 

Pj J dD OV OXj 


So, we infer that Vw is symmetric, divw is constant, and there is a function ip such that ()2.5I) 
holds. 


2.2 Proof of (ii) 

The transmission conditions (|2.81) on dD can be rephrased as 

w|+ = w|_, (j s Vw|_|_^ = <7 c Vw|_ia (2.32) 

Let ti and t 2 be two orthonormal tangent vector fields to dD. Then, we have 
(div w)_ = ((Vw)-v, v) + ((Vw)_ti,ti) + ((Vw)_t 2 ,t 2 ), 


7 






















and 


(div w)+ = ((Vw)|i/, v) + ((Vw)+ti,ti) + ((Vw)+t 2 ,t 2 ). 

Here (divw)_ denotes the limit of divw to dD from inside D , and (divw) + denotes that from 
outside D. Since 

((Vw)-tj,tj) = ((Vw) + tj,tj), j = 1,2, 

we have 

(divw)_ — (divw) + = ((Vw)_i/,i/) — ((Vw) + v,i/). 

It then follows from the second identity in (12.3211 that 


|((Vw) T --(Vw)_)i/,^ = (divw)_ — (divw)+. 

On the other hand, since (Vw) + is symmetric, we obtain 

/((Vw)I - —(Vw)_)i/,tj\ = (*/,(Vw)_tj) - —((Vw )-v,tj) 

\ O s / O s 

= (Vw)+t j) - ((Vw)+i/,t j) = 0. 

We then infer from (I2.33P and (|2.34l) that 

(Vw) T --(Vw)_ ) v = (divw)_i/ — (divw) + zv. 


-x, xeD, 


(2.33) 


Recall that divw is constant in Q\D. Let 

(divw) 


v x = w x — 


2 + £c 
cr s 


(2.34) 


(2.35) 


(2.36) 


Then one can see from (12.3511 that 


Or 


o. 


(Vv) J -(Vv) v — (divv)z^ = 0 on dD. 


(2.37) 


Let g be a smooth vector field on D. It follows from (12.3711 and the divergence theorem that 


= / v ■ f(Vv)g - — (Vv) T g - (divv)g) da 
J dD \ O s J 

= J div ^(Vv)g — — (Vv) T g — (divv)g^j dx. 


One can easily show that 


div (Vv)g- : (Vv) J g — (divv)g = Vv J : Vg--Vv : Vg — (divv)(divg), 

Os J Os 


and so we obtain 


/ Vv J : Vg — — Vv : Vg — (div v) (div g) = 0. 
JD o s 


(2.38) 


Vv := -(Vv + Vv T ) 


and 


Vv := -(Vv — Vv T ) 


Using notation 













it can be rewritten as 


(1 — —) [ Vv : Vg — (1 + —) [ Vv : Vg — [ (divv)(divg) = 0. (2.39) 

a s Jd a s Jd Jd 

If a c > a s , then we take g = v so that 

(1 - —) [ | Vv| 2 - (1 + —) [ | Vv| 2 - f (divv) 2 = 0. (2.40) 

a s jd a s Jd Jd 

Thus, we infer that v is constant in D and hence 


, . (divw) + 

w(x) = --—x + a constant vector, x £ D. 

\ / 2 _|_ <7 c ' 


(2.41) 


If a c = oo, then u is constant on dD , and hence (Vw)t = 0 on dD for any tangential vector t 
to dD. Since Vw is symmetric and divw is constant, it implies that 

(Vw)z' = cqu on dD 

for some constant cq. So, we can see that (ii) holds. 


3 Existence of solutions on confocal ellipsoids 


We first mention that the solution w to m is unique in the sense that if w\ and W 2 are two 
solutions (with different k, A’s, and d’s), then w\ = Cw 2 + E for some constants C and E. 
In fact, if Wj is a solution to (11.41) with k = kj ^ 0, A = A j and d = dj (j = 1,2), then 
vj = w 1 — satisfies Aw = 0 in O \ D and Vw = 0 on <917, so we have that w must be a 


constant. 

We now construct a solution to (jl.4|) when D and 12 are confocal ellipsoids. To do so, assume 
that dD is given by 


2 2 2 

X 1 . x 2 1 x 3 _ i 

„2 ' „2 ' „2 

C 1 c 2 c 3 


We then use the confocal ellipsoidal coordinates p, p, £ such that 


(3.1) 


c\ + p 


+ 


4 + p 


+ 


ci + p 


J2 2 

Ju 1 Ju O 

1 + + 


c{ + p C 2 + p C 3 + p 


t 2 t 2 

Jb 1 Jbrx 

+ + 


cf + e 






1, 

1, 

1, 


subject to the conditions — c 2 < £ < — c| < p < —c 2 < p. Then the confocal ellipsoid <912 is 
given by p = po for some p 0 > 0. 

Let 

dip) = ( c i + p)(, c 2 + P)( c 3 + P)i (3-2) 

and dehne 


<Pj(p) = 


1 


'p (c? + s)y/g(sj 


ds, j = 1,2,3. 


(3.3) 
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Then the function w defined by 


w(x) = 


1 1 1 3 1 3 

-2 l 7m iS ~2 + 


(3.4) 


is a solution of (11.41) . In fact, we can see that 

3 


_d_ 

dxi 




#) 


5/9 


*5 I &7 - * Wli ' 


Since 


we have 


9=1 


=“ 


V - T j = 1 

“t (c? + p)\/5(p) V 7 ®’ 


_d_ 

dxi 




= -ipi(p)Xi, 


from which we see that 

Vw(x) = -(<Pl(p)xi,<p 2 (p)x2,(p3(p)X3) + (V1(P0)X1,IP2(P0)X2,<P3(P0)X3)- 
Using the relation 


dp 2 Xi 


+ 


Xn 


+ 


-1 


(3.5) 


(3.6) 


d Xl cj + p[{c\ + pY {cl + pY (4 + pY\ ’ 

we see that Aw is constant. Note that \7w = 0 on 5fl (p = p o) and \7w = Ax on dD where 

A = diag[<£>i(p 0 ) - <^i(0),^2(Po) - <£ 2 ( 0 ), <Pz{po) - ^3(0)]- (3-7) 

We emphasize that A is negative-definite. 


4 Proof of Theorem 11.11 

Let w be the solution to (|1.4j) with A = cl. We notice that c / 0. Indeed, if c = 0, then we have 


0^fc|ll\D|= / A w dx = 

Jn\D Jan 


f (L d(7 — [ ^<4t = 0- / „.d<4r = 0, 
Jan dv Jqd du J dD 


which is a contradiction. Since c ^ 0, by introducing new variables 


y = x + -d, 

c 


we may assume that d = 0. Set 


„ d d 

Aij — Xj — Xi - 

dxi dx 


l n I * Y 3- 
i OXj 


(4.1) 


It is worth mentioning that A tJ is the angular derivative. Observe that A^ commutes with A, 
namely, A^ A = A Aij. So, we have A AijW = 0 in 0 \ D. Note that AijW = 0 on dQ. Since 
Vtc(x) = cx on dD, we see that A^jw = 0 on dD. Then the maximum principle yields that 


A^w = 0 in O \ D. 


(4.2) 
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Since Aw = k in 12 \ D, w satisfies the ordinary differential equation 

d 2 w 2 dw . 

TT + - 7 - = i: m 12 \ 14 
or z r or 

for r = |x|. Choose a ball B with B C 14 \ D. By 114.31) . w is of the form 


k 


w(r ) = —r 2 H-(- /c 2 in B 

6 r 

for some real constants k\ and k 2 ■ Since 14 \ D is connected and 

k 2 k\ 


we have from (14.41) 


A ui-r 2 - &2 ) = 0 in 14 \ D, 

'or' 


u;(r) = —r 2 + — + k 2 in 14 \ 14. 
6 r 


Since ^ = 0 on 914, we must have 


k kl n on 

—r it = (J on ail, 


(4.3) 


(4.4) 


(4.5) 


and hence 


3 

" = —— on 912. 


This means that 914 = 8Br{ 0) for some R > 0. Therefore we have 

_ , . k kR? x . — 

Vw(x) = -x-—- 3 , x € 14 \ 14. 

Since Vrc(x) = cx for all x E 914, we must have 

/c fci? 3 1 

3 3~7 3 

or r = constant for all x E 9Z4. It means that 9Z4 is a sphere centered at 0 . This completes the 
proof. □ 


= c on 914, 


5 Newtonian potential formulation 


In this section we reformulate the problem (11.41) in terms of the Newtonian potentials and relate 
the problem with known characterization of ellipsoids using the property of the Newtonian 
potential. 

Suppose that (11.41) admits a solution and let w be the solution. Notice that by the second 
equation of (11.41) w is constant on each connected component of 914, and by the third equation 
of (11.41) u;(x) = |x • Ax + d • x + C for x E 9Z4 for some constant C. Fix x ^ 14 \ 14. We obtain 
from the divergence theorem that 


k[ r(x — y)dy = f [Aiu(y)T(x - y) - w(y)A y r(x - y)] dy 
Jn\D Jn\D 


L 

L 


IQ\D 

A(y)r(x-y)-»(y)Ar(x-y) 

{V • Ay + 1 / • d)T(x - y ) _ (ly . Ay + d • y + C)^-T( X - y) 


dcr (y) ~ [ w(y)^-T(x-y)da(y) 
JdQ vvy 

dcr(y) 


9za, 


ion 


9 

w( y ) 7 ^ r ( x - y)dcr{y). 
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If x € M 3 \ fl, then by dealing with the last integral for each component of dfl we have 

[ w (y)-J— T (x-y) dcr (y) =°> 

Jan av y 

and hence 

k [ T(x - y)dy = -Tr A f T(x - y)dy. 

Jn\D Jd 

We can find a relation between k and Tr A from this formula. In fact, we have 

lim [ r(x-y)dy = |fi\£>| 

Ix| yoo I (xj Jn\D 

and a similar formula for J D T(x — y)dy. So we obtain 

k\n\D\ = -TrA|D|. (5.1) 


We then see that 


k / T(x — y)dy + Tr A / T(x - y)dy = k\Q\ Afo (x) - IVd( x) 
Jn\D Jd l 


(5.2) 


where A^ and Ad are the (averaged) Newtonian potentials on £1 and D, respectively, namely, 


Sn(x):= MX r(x ~ y)<iy ’ 


(5.3) 


and similarly for Ad- 

If x E D, then we have for some constant C* 


k f T(x - y)dy = -Tr A [ T(x - y)dy + ^x • Ax + d • x + C* 
Jfi\D Jd 2 


and hence 


1 


k\Q\ Aq(x) - Ad(x 
We have shown that if (11.41) admits a solution, then 


= —x • Ax + d • x + C*. 
2 


An(x) - Ad(x) = 


o, x e m 3 \ n, 

a quadratic polynomial, x £ D. 


(5.4) 


(5.5) 


So we may reformulate the question: If (|5.5j) holds, then D and Q, are confocal ellipsoids. This 
is reminiscent of a question related to the Newton potential problem: If a Newtonian potential 
of a simply connected domain is a quadratic polynomial in the domain, then the domain must 
be an ellipsoid. This problem has been solved by Dive [5] and Nikliborc [16] (see also |3| and 

m)- 
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